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(57) ABSTRACT 

Methods and systems for converting an image contrast evo- 
lution of an object to a temperature contrast evolution and 
vice versa are disclosed, including methods for assessing an 
emissivity of the object; calculating an afterglow heat flux 
evolution; calculating a measurement region of interest tem- 
perature change; calculating a reference region of interest 
temperature change; calculating a reflection temperature 
change; calculating the image contrast evolution or the tem- 
perature contrast evolution; and converting the image contrast 
evolution to the temperature contrast evolution or vice versa, 
respectively. 
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METHODS AND SYSTEMS FOR 
MEASUREMENT AND ESTIMATION OF 
NORMALIZED CONTRAST IN INFRARED 
THERMOGRAPHY 

5 

I. CROSS REFERENCE TO RELATED 
APPLICATIONS 

This application claims the benefit of the filing date of U.S. 
provisional patent application No. 61/293,426 incorporated 10 
herein by reference, which was filed on Jan. 8, 2010 by the 
same inventor of this application. 

II. ORIGIN OF THE INVENTION 

15 

The invention described herein was made by an employee 
of the United States Government and may be manufactured 
and used by or for the Government of the United States of 
America for governmental purposes without the payment of 
any royalties thereon or therefor. 20 

III. FIELD OF THE INVENTION 

The invention is directed to the fields of non-destructive 
evaluation (NDE) and of processing of data acquired from an 25 
infrared camera. Specifically, the systems and methods 
described herein concern the performance of a series of steps 
for analyzing and processing digital data comprising a plu- 
rality of infra-red (IR) video images acquired by a system 
used for non-destructive evaluation. 30 

IV. BACKGROUND OF THE INVENTION 

The embodiments described herein relate generally to the 
field of non-destructive evaluation (NDE) using thermo- 35 
graphic images. Infrared (IR) flash (or pulsed) thermography 
is an example of a technique for NDE used in the inspection 
of thin materials such as laminated or bonded composites in 
the aerospace industry. IR flash thermography is used to 
detect delamination-like anomalies, although other anoma- 40 
lies, such as surface cracks, may be detected. 

An example of hardware equipment for an IR flash ther- 
mography system comprises a flash lamp (source of light/ 
heat), a flash hood or housing, a flash power supply/trigger 
unit, a flash duration controller, an IR camera for capturing 45 
video images, data acquisition electronics, and a computer. 
The computer may be used for controlling the flash trigger, 
for acquiring video data from the IR camera, for displaying 
data, and for post -processing of the acquired data. 

In one example of an NDE technique using IR flash ther- 50 
mography, a single sided or reflection technique is used 
wherein the flash lamp (heat source) and the IR camera (de- 
tector) are on the same side of a test object undergoing inspec- 
tion. A plate is provided as the test object with a round 
delamination in the center. After applying heat to the top 55 
surface of the test object by triggering the flash lamp, the top 
surface area surrounding the anomaly cools faster than the top 
surface (footprint) area above the anomaly. The IR camera 
captures a sequence of images of the surface temperature in 
terms of pixel intensity and represents the anomaly as a hot 60 
spot (e.g., an area warmer than the surrounding area or the 
reference region of interest (ROI)). The hot spot is about the 
size and shape of the anomaly footprint. Relative pixel inten- 
sity, i.e., the difference in pixel intensity between the hot spot 
(measurement ROI) and the surrounding area (reference 65 
ROI), varies with the post-flash time. Deeper anomalies 
appear in the IR video data at later times compared to the near 
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surface anomalies. After the initial appearance of an anomaly 
in the IR video data, the relative pixel intensity continues to 
increase with time. The relative pixel intensity of the anomaly 
reaches a peak at a certain time, and then the relative pixel 
intensity decays until the temperature of the indication area 
and the temperature of the surrounding area become equal. 

V. SUMMARY 

In one aspect, disclosed is a method for converting an 
image contrast evolution of an object to a temperature con- 
trast evolution, the method comprising: calculating a mea- 
surement region of interest temperature change AT; calculat- 
ing a reference region of interest temperature change AT,.^, 
calculating a reflection temperature change AT re ^; calculat- 
ing the image contrast evolution C„!\ and converting, using a 
processor, the image contrast evolution to the temperature 
contrast evolution. 

In another aspect, disclosed is a method for converting a 
temperature contrast evolution of an object to an image con- 
trast evolution, the method comprising: calculating a mea- 
surement region of interest temperature change AT; calculat- 
ing a reference region of interest temperature change AT re/ ; 
calculating a reflection temperature change AT^; calculat- 
ing the temperature contrast evolution U; and converting the 
temperature contrast evolution to the image contrast evolu- 
tion. 

In another aspect, disclosed is an apparatus for converting 
an image contrast evolution of an object to a temperature 
contrast evolution, the apparatus comprising one or more 
processors and one or more memory units coupled to the 
processors. The apparatus is configured and arranged to cal- 
culate a measurement region of interest temperature change 
AT, to calculate a reference region of interest temperature 
change AT re;/ , to calculate a reflection temperature change 
AT ;iy?) to calculate the image contrast evolution C^, and to 
convert the image contrast evolution to the temperature con- 
trast evolution. 

In another aspect, disclosed is an apparatus for converting 
a temperature contrast evolution of an object to an image 
contrast evolution, the apparatus comprising one or more 
processors and one or more memory units coupled to the 
processors. The apparatus is configured to: calculate a mea- 
surement region of interest temperature change AT; calculate 
a reference region of interest temperature change AT re y; cal- 
culate a reflection temperature change AT^; calculate the 
temperature contrast evolution C ! \ and convert the tempera- 
ture contrast evolution to the image contrast evolution. 

In another aspect, disclosed is a method for assessing emis- 
sivity of an object being inspected by an infrared flash ther- 
mography system, the method comprising the steps of: mea- 
suring a pre-flash temperature at a measurement region of 
interest W°; calculating a camera constant C cam ; and calcu- 
lating the emissivity of the object according to the equation 
of: 


In yet another aspect, disclosed is a method for calculating 
an afterglow heat flux evolution S post fl ash , the method com- 
prising the steps of measuring a reflection temperature T re ^; 
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and calculating the afterglow heat flux evolution S> posT ft as/ , 
according to the equation of: 

S^efart^orCT^Z-T^ 04 ), 

wherein T re ^° is the reflection temperature at time of flash and 
ct is the Stefan-Boltzmann Constant. 

VI. BRIEF DESCRIPTION OF THE DRAWINGS 

to 

FIG. 1 is a schematic diagram of an exemplary single-sided 
flash thermography system in accordance with some embodi- 
ments described herein. 

FIGS. 2(a) and ( b ) are a schematic diagram of a plate 
depicting a gapping delamination with measurement and ref- 1 5 
erence regions of interest (ROIs) and a single infrared (IR) 
video image of a test object with back drilled flat bottom 
holes, respectively, in accordance with some embodiments 
described herein. 

FIG. 3 is an infrared image of a plate with back drilled flat 20 
bottom slots along with foils and tape (post-flash time= 0.6 
sec), in accordance with some embodiments described 
herein. 

FIG. 4 is a graph of the value of pixel intensity as a function 
of position along a horizontal cross line in a frame taken at 0.6 25 
seconds after the flash, in accordance with some embodi- 
ments described herein. 

FIG. 5 is a graph of the pixel intensity evolution as a 
function of time on the foil ROI (lower evolution from the left 
foil ROI and upper evolution from the right foil ROI) with the 30 
peak time of 0.6 seconds indicated by the dashed vertical line, 
in accordance with some embodiments described herein. 

FIG. 6 is a graph of the pixel intensity evolution as a 
function of time for the tape ROI, in accordance with some 
embodiments described herein. 35 

FIG. 7 is a graph of the ROI intensity evolutions for the 
measurement (upper evolution) ROI and the reference (lower 
evolution) ROI, in accordance with some embodiments 
described herein. 

FIG. 8 is a graph of simple contrast evolution at the mea- 40 
surement ROI with the peak time of 0.5 seconds indicated by 
the dashed vertical line, in accordance with some embodi- 
ments described herein. 

FIG. 9 is a graph of computed evolution of the reflection 
temperature, in accordance with some embodiments 45 
described herein. 

FIG. 10 is a graph of estimated temperature evolutions for 
the measurement ROI (middle curve), the reference ROI (bot- 
tom curve), and the reflection or the background ROI (top 
curve), in accordance with some embodiments described 50 
herein. 

FIG. 11(a) is a graph of normalized image contrast com- 
puted using the pixel intensity (upper curve) and one fourth 
power of the pixel intensity, in accordance with some embodi- 
ments described herein. 55 

FIG. 11 (b) is a graph of measured normalized image con- 
trast (lower curve) and estimated normalized temperature 
contrast (upper curve), in accordance with embodiments 
described herein. 

FIG. 12 is an exemplary screen shot displaying normalized 60 
temperature contrast evolution (upper curve) and normalized 
image contrast evolution (lower curve) computed using soft- 
ware that implements simulation and data analysis methods 
described herein in accordance with some embodiments. 

FIG. 13 is a graph of the pixel intensity to temperature 65 
contrast ratio, in accordance with some embodiments 
described herein. 
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FIG. 14 is a graph of simulated temperature evolutions with 
an assumed reflection temperature evolution, in accordance 
with some embodiments described herein. 

FIG. 15 is a graph of estimated image contrast and of 
simulated temperature contrast, in accordance with some 
embodiments described herein. 

FIG. 16 is a graph of estimated afterglow flux as a function 
of time charted from the first post-flash frame, in accordance 
with some embodiments described herein. 

FIG. 17 is a graph of flash flux evolution as a function of 
time until the first post-flash frame, in accordance with some 
embodiments described herein. 

FIG. 18 is a flow diagram illustrating methods for convert- 
ing an image contrast evolution of an object to a temperature 
contrast evolution, for modeling afterglow flux, and for 
assessing the test object emissivity, in accordance with some 
embodiments described herein. 

FIG. 19 is a flow diagram illustrating a method for convert- 
ing a temperature contrast evolution of an object to an image 
contrast evolution, in accordance with some embodiments. 

Because the embodiments described herein are subject to 
various modifications and alternative forms, it should be 
understood that the drawings and detailed description are not 
intended to limit the invention to the particular embodiments 
described herein. Instead, the following description is 
intended to cover all modifications, alternatives, and equiva- 
lents falling within the scope of the invention as defined by the 
appended claims. 

VII. DETAILED DESCRIPTION 

The objects in the drawings are not necessarily to scale and 
certain features may be shown exaggerated in scale or in 
somewhat generalized or schematic form in the interest of 
clarity and conciseness. In the description which follows like 
parts may be identified throughout the specification and 
drawings with the same reference numerals. The foregoing 
brief description of the drawings is provided for a more com- 
plete understanding thereof. It should be understood, how- 
ever, that the embodiments described herein are not limited to 
the precise arrangements and configurations shown. 
Although the design and use of various embodiments are 
discussed in detail below, it should be appreciated that the 
embodiments described represent concepts that may be 
embodied in a wide variety of contexts. The specific aspects 
and embodiments described herein are merely illustrative, for 
it would be impossible or impractical to include all of the 
possible embodiments and contexts in this description of the 
invention. Upon reading this description, alternative embodi- 
ments will be apparent to persons of ordinary skill in the art. 

During the processing of digital data acquired from an 
infrared camera, normalized image (or pixel intensity) con- 
trast and normalized temperature contrast differ in value for 
objects having a value of emissivity other than one. There- 
fore, for more accurate processing of data obtained from 
objects having an emissivity value other than one, the two 
quantities should not be treated as the same. In the co-pending 
U.S. patent application Ser. No. 12/900,644, filed Oct. 8, 
2010 and entitled “Methods and Systems for Characterization 
of an Anomaly Using Infrared Flash Thermography”, which 
is hereby incorporated in its entirety by reference, Koshti 
distinguishes and defines the two quantities of normalized 
pixel intensity contrast and normalized temperature contrast. 
Embodiments described herein establish formulas to estimate 
the normalized temperature contrast from the normalized 
pixel intensity or image contrast using the methods and sys- 
tems described herein. The method includes steps for the 
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measurement of the reflection temperature evolution by com- 
paring the simulation temperature contrast with the measured 
pixel intensity contrast. The method also comprises the act of 
estimating the incident heat flux. Ideally, the simulation 
should model the compound heat source flux evolution which 
also includes the thermal afterglow. The effect of the reflec- 
tion temperature in the pixel intensity should also be 
accounted for to seek a more refined estimation of the tem- 
perature contrast profile from the pixel intensity evolution 
data. 

The method of measuring reflection temperature evolution 
is described herein. F irst, the reflection temperature evolution 
is established based on infrared (IR) flash thermography (or 
simply any transient thermography) data acquisition. This 
acquisition of data requires a test object, a high emissivity 
tape with known emissivity, and a diffused but highly reflec- 
tive metal foil with known reflectivity. The method also 
records the steady state pre-flash temperature of the object 
using a they thermocouple (or other contact sensor) or an 
accurate radiometer. The IR datacube is recorded after per- 
forming the technique of IR flash thermography. 

The method further comprises data acquisition from four 
regions of interests (ROI). One region is for the measurement 
area. The second ROI is for the reference area. The third ROI 
is for the measurement of the reflection temperature. The 
fourth ROI is for the measurement of the pre-flash tempera- 
ture of the high emissivity tape. 

Using formulas given herein, the method estimates the 
reflection temperature evolution. Then, the method computes 
the temperature contrast from the IR data. The emissivity 
factor is defined to relate the temperature contrast to the 
image contrast. 

The method uses the reflection temperature evolution to 
model the afterglow flux of the flash source. Using the esti- 
mated compound source evolution in the simulation software, 
the method can estimate the temperature contrast evolutions 
and then estimate the image contrast profiles on simulated 
voids. 

The invention provides an emissivity estimation technique 
using the same IR camera. The technique estimates the emis- 
sivity of an object for a desired thermal wavelength. It is 
shown that using the foil-tape (or tape-foil) technique during 
the IR shot, the transient reflection temperature or the reflec- 
tion temperature evolution can be recorded. If the IR camera 
is programmed to use the real-time reflection temperature 
using the formulas described herein, it can provide real-time 
object surface temperature. The IR camera can be pro- 
grammed to estimate the object emissivity in real-time using 
the formulas derived herein for the foil-tape technique. 

Referring now to the drawings, FIG. 1 is a schematic dia- 
gram of an exemplary single-sided flash thermography sys- 
tem in accordance with some embodiments described herein. 
As shown in FIG. 1, the equipment for an infrared (IR) flash 
thermography (IRFT) system 105 in accordance with at least 
one of the exemplary embodiments described herein com- 
prises a flash lamp (source of light/heat) 110, a flash-hood 
115, at least one flash power supply/trigger unit 120, at least 
one flash duration controller 125, an IR camera 130 for cap- 
hiring video images, camera data acquisition electronics 135, 
and a computer 140. The computer 140 is used for controlling 
the flash trigger unit 120 forpowering the flash lamp 110, data 
acquisition of the camera video data, data display, and post 
processing of the acquired data. Hie flash-hood 115 may be 
made from sheet metal and have highly reflective (e=0.05) 
inside wall surfaces. The flash-hood may form a box-like 
structure, although other structures that form an interior vol- 
ume, such as a cylinder, hemi-sphere, etc. may also be used. 
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However the flash-hood is structurally configured, one por- 
tion of the hood has a laige opening so as to enable the hood 
to be positioned over a test obj ect 1 45 . The side opposite to the 
opening has a hole in the center to provide a window for the 
5 lens of the IR camera 130 that is positioned outside of the 
flash-hood 115. The IR camera 130 is focused at a surface of 
the test object 145 positioned within the hood. At least one 
flash lamp 110 (two are shown in FIG. 1) is located within the 
inner wall of the flash-hood 115 in proximity to the IR camera 
to 130. The flash lamp(s) 110 is covered with a glass shield 
which is highly transmissive for light but highly emissive for 
thermal radiation. The flash lamp(s) 110 when powered 
directs illumination towards the test object 145 without 
directly shining light upon or into the IR camera 130. The 
15 flash-hood 115 functions as a housing that contains most of 
the intense flash from the flash lamp 110. 

If the size of the test object 145 can be accommodated 
inside the hood 115, then the test object 145 is located at the 
hood opening or slightly inside the hood 115. Otherwise, the 
20 test object 145 is located slightly outside of the hood 115 
opening. A short duration (e.g. 3 millisecond), intense (12 kJ) 
flash is triggered by the computer 140. Data acquisition is 
initiated a few seconds before the flash, and it continues until 
a prescribed duration of time has expired. The camera 130 
25 provides a sequence of IR images (or frames) called the 
data-cube (or digital video) of the test object surface taken at 
the chosen frame rate (e.g. 60 Hz or 60 frames per see). As 
described in further detail herein, the intensity (numerical 
value) of each pixel of the image is a function of the surface 
30 temperature of the corresponding area observed of the test 
object 145 at the time of the image frame. The flash causes the 
surface of the test object 145 to warm up slightly, and the heat 
then starts to dissipate rapidly. The surface cools through 
thermal radiation, convection, and conduction. Heat conduc- 
35 tion within the test object 145 is considered to be the domi- 
nant heat transfer mode until the temperature gradients within 
the test object become small. At later times, the heat conduc- 
tion is of the order of the combined effect of heat convection 
and radiation. The IR data acquisition and data analysis uti- 
40 lizes the thermal data in the short duration immediately after 
the flash where the thermal dissipation is dominated by the 
heat conduction within the test object 145. 

Heat exchange across the boundaries due to convection can 
be assumed to be zero if the Biot number (N Si =hL/k) is less 
45 than 0. 1 . For example, a half centimeter thick graphite/ epoxy 
(k=0.64 W/mK) plate, using h= 1 0 W/m 2 K, has a Biot number 
of 0.078. Therefore, heat conduction is the dominant mode of 
heat transfer in this graphite/epoxy plate example. Thinner 
test objects tend to equalize the temperature within the test 
50 object very quickly and have a relatively longer cooling time 
by heat loss to environment. 

FIG. 2(a) is a schematic diagram of a plate serving as the 
test object 145. FIG. 2(a) depicts a gapping delamination or 
anomaly 220 with measurement and reference regions of 
55 interest (ROIs) in accordance with embodiments described 
herein. FIG. 2(b) shows a single infrared (IR) video image of 
a test object with a flat bottom hole in accordance with at least 
one embodiment. Referring now to FIGS. 1 and 2(a), the test 
object 145 is made of a thermally isotropic material with 
60 constant thickness 205 that fits inside the hood 115. The test 
object 145 is supported at the corners on insulating standoffs 
150 and the hood is oriented vertically thereabove. Assuming 
that the flash intensity is uniform over the test object top 
surface 210, the heat conduction will be in a direction normal 
65 to the test object top surface 210 in most of the acreage area 
(i .e., area away from edges of the test obj ect and flash bound- 
ary )• The heat is conducted uniformly from the top surface of 
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the test object 210 to the bottom surface of the test object 215. 
Normal heat conduction will be obstructed by an anomaly 
such as a small round gapping delamination 220 at the center 
of the test object, as depicted in FIG. 2(a). The delamination 
220 is at some depth (d) 221 below the top surface of the test 
object 210. The delamination may also be characterized as 
having some diameter (D) or width 222 and a gap 223. The 
volume hounded by the anomaly 220 on one side and the top 
surface 210 on the other side is called the heat trapping 
volume 225. 

The top surface area 230 surrounding the anomaly 220 
cools faster than the top surface (footprint) area 235 above the 
anomaly 220. Hie IR camera 130 captures the surface tem- 
perature image in terms of the pixel intensity and shows the 
anomaly 220 as a hot spot (e.g. an area warmer than the 
surrounding area) which is about the size and shape of the 
anomaly footprint. The relative pixel intensity of the hot spot 
changes with the time. Deeper anomalies appear at later times 
in the IR video compared to the near surface anomalies. After 
the appearance of an anomaly in the IR video, its relative pixel 
intensity continues to increase with time. The relative pixel 
intensity associated with the anomaly reaches a peak at a 
certain time, and then the relative intensity decays until the 
indication area temperature and the surrounding area tem- 
perature become equal. 
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-continued 



where T Ju;/ is the top surface temperature, T amb is the ambient 
temperature of the air contacting the surface, q, ; is the heat 
convection rate (cal-sec -1 or BTU-hr -1 ^ h com is the heat 
convection film coefficient (cal-sec -1 - 0 C. -1 -cm -2 , BTU-hr -1 
F - 1 -ft -2 or W-K - 1 -m -2 ), A is the surface area of the plate(cm 2 , 
m 2 or ft 2 ), and R ; , is the convective resistance (sec-cal -1 - 0 C. 
orhr-BTU -1 -°F.). Fornatural (unforced) convection, atypical 
value for the heat convection film coefficient is 4 W-K - 1 -m -2 . 
Using equations (3)and (4), the heat convection rate across 
each side of a delamination air gap with one surface at tem- 
perature Tj and other surface at temperature T 2 (r i>T 2 ) is 
given by equation (5a), 

Tl - T am b Tamb - T; (5a) 

qh allies A) - (1 !K om A) 

Further, the ambient temperature of the air for thin gaps is 
given by equation (5b), 


FIG. 2(b) shows one of the frames in a sample IR data-cube 
from a reinforced carbon-carbon (RCC) part with round flat 
bottom holes (FBFls) machined into the back surface to simu- 30 
late gapping delaminations of differing depths. IR indications 
240 from different fiat bottom holes are clearly visible in the 
IR image. A measurement ROI 245 is shown as a square in the 
middle of one of the anomalies. A reference ROI 250 is shown 
as a square outside the anomaly area. As the test concludes, 35 
the test object continues to cool down to ambient temperature 
through heat convection and radiation. 

In some embodiments, such as for a test object that is a flat 
plate with a large diameter, the top surface may be represented 
as being at a temperature T : and similarly the bottom surface 40 
is at a temperature T,. A heat conduction rate for an isotropic 
slab may accordingly be given by equations (1 ) and (2), 
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Ti-Tz 
C L/kA ) 


and 



(1) 45 

(2) 

50 


where q* is the heat conduction rate (cal-sec -1 or BTU-hr -1 ), 
k is the thermal conductivity (cal-sec -1 - 0 C. -1 -cm -1 or BTU- 
lir -1 - 0 F. -1 -ft -1 ), A is the cross sectional area of the plate (cm 2 
or ft"), L is the thickness of the plate (cm or ft), and R A . is the 
conductive resistance (sec-cal -1 - 0 C. orhr-BTU -1 -°F.). Equa- 
tions (1 ) and (2) may be used for representing heat conduction 
across a delamination gap L with one surface at temperature 
Tj and the other at temperature T 2 . 


Ti + T 2 (5b) 


When the air gap becomes smaller than a threshold value 
represented as the Biot number, the amount of heat exchange 
by convection is replaced by heat conduction. The heat 
exchange across the boundaries due to the convection can be 
assumed to be zero if the Biot number is less than or equal to 
0.1 (N Sl =hL/ks0.1 ). This relationship gives an expression for 
the threshold air gap thickness as L,,, r =0. 1 k/h. Taking the air 
conductivity to be 0.026 W/mK and assuming the convection 
coefficient to be 5 W/m 2 K, then a threshold gap thickness of 
0.020 inch (0.5 nun) is calculated, beyond which heat con- 
vection is applicable. Heat conduction is present if the gap 
thickness is less than the threshold gap thickness. The nar- 
rowing of the air gap results in increased heat conduction 
across the air gap, which manifests as a decrease in the peak 
contrast value of the void. Loss of the peak contrast value by 
20% may not occur until the gap narrows significantly com- 
pared to the threshold value of the gap thickness. 

In some embodiments, such as for a flat plate with a large 
diameter, the radiative heat transfer rate is given by equations 
(6), (7), and (8), 


T sur f ~ T re ji 

9r ~ (l/M) ’ 

h r = ea-(T surf + T„ fl )(T 2 urf + T 2 rljl ), and 


( 6 ) 

(7) 

( 8 ) 


In some embodiments, such as for a flat plate with a large 60 
diameter, the heat convection rate from a surface is given by 
equations (3) and (4), 



where T^ is the apparent temperature experienced by the 
object surface due to the camera-side background tempera- 
hire in units of Kelvin or Rankine, T amb is the temperature of 
the air surrounding the object in Kelvin or Rankine, cr is the 
Stefan-Boltzmann Constant 5.670400xl0 -8 Wm -2 K -4 , e is 
the emissivity of the test object surface (assuming same emis- 
sivity on top (front) and bottom (rear) surfaces), h r is the heat 
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radiation transfer coefficient (cal-sec _1 -° C. -1 -cm -2 or BTU- 

hr _1 -° F. -1 -fT 2 ), and R, is the radiative resistance (sec-cal _1 -° 

C. or hr-BTU _1 -° F.). At a temperature of 300 K for a 10 K 
temperature difference between the object surface tempera- 
ture and the corresponding reflection (background) tempera- 5 
ture, the radiation transfer coefficient is less than 6 W-K -1 - 
m -2 . 

In some embodiments, such as for a large but thin slab of 
isotropic material with uniform thickness that is exposed to a 
heat pulse on the front surface, the resulting one dimensional to 
transient heat transfer is described by equations (9) and (10a), 


dT d 2 T 

(9) 

- dl =a 'd ? * nd 

15 

II 

o 

II 

S' 

(10a) 


where a is the test object (slab) diffusivity (cm 2 /s), and z is the 
distance along the thickness direction. The total heat transfer 20 
rate at the slab surface is given by equation (10b), 

q=q*+q*+qr- (iob) 

Neglecting the air emission, the heat flux boundary condi- 
tions are given by equations (11) and (1 2), 25 



? — h f ron t (T front 


: 7 hrear(T r ear 


' Tomb— front ) ^(Tfront ^refl- front) 


Tomb-rear ) ^ £ )Trear Trcfl-rair)' 


01) 

( 12 ) 


30 


where T(t=0)=T,„ is the initial steady state temperature of the 
test object before the flash and h fronr and h rear are the convec- 35 
tion coefficients at the front and rear surfaces, respectively. 

In practice, an IR camera measures pixel intensity. Radio- 
metric cameras can measure the surface temperature. Radio- 
metric cameras are not known to be set-up for compensating 
the transient surrounding temperature. Radiometric cameras 40 
may have input for a constant reflection (background) tem- 
perature. They may use a reflective foil to assess the stable 
reflection temperature. For the exemplary embodiments 
described herein, either a radiometric IR camera (i.e. measure 
apparent temperature) or an IR camera that measures the 45 
irradiance (i.e. measure pixel intensity) may be used because 
the exemplary embodiments use a method described herein to 
capture the transient reflection temperature. The method 
assumes that camera pixel intensity response is proportional 
to the camera irradiance in the desired temperature range. 50 

Considering the definition of the normalized contrast based 
on the pixel intensity, its relationship with the temperature of 
the test object will be derived. The pixel refers to the camera 
image picture element. The pixel size refers to the instanta- 
neous field of view (IFOV), which is the corresponding area 55 
of the object surface imaged in the pixel. A selected pixel grid 
area, usually in the form of a rectangle, is called the image 
region of interest (ROI, or when plural ROIs). The size of the 
ROI then refers to the area of the object surface imaged in the 
ROI. The corresponding area of the test object is referred to as 60 
the object ROI. The camera measures total radiation incident 
on its detector array element and displays it as a pixel in the 
image. The detector element response or the pixel intensity 
registered by the camera is governed by the following equa- 
tion (13a), 65 


(13a) 


10 

where W is the average pixel intensity due to the heat irradi- 
ance (measured in units of W/m 2 , cal-sec _1 -cm -2 , or BTU- 
hr _1 -fT 2 ) of the corresponding camera detector elements 
from the object ROI measured in gray scale bit value (positive 
integers). The right hand terms in equation (13a) are given by 
equations (14), (15), and (16), 


W oij .=€TW' oij , 

(14) 

W r ^=(l-e)rW' ' re jj, and 

(15) 


(16) 


where Vf obJ is the contribution to the pixel intensity due to the 
heat emission from the object ROI, is the contribution to 

the pixel intensity due to the heat reflection from the object 
ROI, and W aw is the contribution to the pixel intensity due to 
the heat emission from the air between the object ROI and 
camera. If the method assumes a perfect focus and that the 
ROI has a uniform temperature at any time, then equations 
(17), (18), and (19) will apply, 


W'^=C„ 

™oT*, 

(17) 

w ' reJ fC ca 

and 

(18) 



(19) 


where W' obJ is the uncompensated contribution to the pixel 
intensity due to the heat emission from the test object ROI, 
W' ref) is the uncompensated contribution to the pixel intensity 
due to the heat reflection from the test object ROI, is the 

uncompensated contribution to the pixel intensity due to the 
heat emission from the air between the camera and object, T 
is the surface temperature in Kelvin or Rankine, T refl is the 
apparent temperature experienced by the object surface due to 
the camera-side background temperature in Kelvin or Rank- 
ine, T ann i s the temperature of the air between the camera and 
theobject in Kelvin or Rankine, C cam aisthecarnera constant, 
€ is the themial emissivity of the test object surface, and x is 
the themial radiation (at camera detection wavelength) trans- 
missivity of the air between the test object and the camera. 
The thermal emissivity is dependent on the thermal wave- 
length for a non-gray body. The wavelength spectmm is 
affected by the surface temperature. Assuming that the emis- 
sivity is high and the test object is a gray body, the emissivity 
is not dependent on the wavelength. However, the IR camera 
measures radiation in a selected band of the wavelength (e.g. 
3-5 micron). Equations (17), (18), and (19) apply to the total 
radiation integrated over all wavelengths. Thus, in reality, in 
equations (17), (18), and (19), the exponent (or power) ofthe 
temperature variable T may not be equal to four. The exponent 
of the temperature variable may be between 2 to 17. The 
following analysis is applicable for these values of the expo- 
nent. 

A calibrated camera provides a linear response with the 
thermal radiation incident on the detector array element. This 
relationship implies that, in a hypothetical case of absolute 
zero temperature, if the camera calibration is extrapolated, all 
detector elements shall provide zero pixel intensity. The 
method assumes that the camera constant is the same within 
a range of 280 K to 350 K. The pixel intensity or the pixel 
grayscale value also depends upon the pixel bit resolution 
(typically 14 bit). The pixel intensity has an upper saturation 
limit dictated by the pixel bit resolution. A typical InSb mid- 
wave camera used in flash thermography operates at a wave- 
length of about 3 to about 5 micron. The transmissivity of air 
is a function of the distance, humidity, air composition, and 
wavelength. A pronounced dip in the transmissivity occurs at 
a wavelength of 4.2 micron due to CO, absorption. Due to a 


w-w^+w^+w^ 
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short distance (~1 ft) between the camera and the test object, 
the contribution of air emission (emissivity~0) is very small 
and is neglected herein. Correspondingly, the air transmission 
(transmissivity~l) is almost 100%. The microbolometer 
long -wave camera works in the 9 to 1 2 micron wavelength and 
has better air transmissivity even for longer distances of sev- 
eral feet (e.g. 30 ft). Therefore, equation (13a) can be simpli- 
fied to equation (13b), 

W«W^.+W^. (13b) 

The normalized image contrast is defined based on the 
pixel intensity as expressed in equations (20), (21), and (22), 


, AW-AW^ 
w AW+AW„ f ’ 

(20) 

A W = W-W°, and 

(21) 

AW r , f = W„ f - W° f , 

(22) 


where is the normalized IR image (pixel intensity) con- 
trast, AW is the change in the pixel intensity of the measure- 
ment ROI after the flash, W is the pixel intensity at the mea- 
surement ROI at the post-flash time t, W° is the pixel intensity 
at the measurement ROI before the flash, AW, c/ is the change 
in the pixel intensity of the reference ROI after the flash, 
is the pixel intensity at the reference ROI at the post-flash time 
t, and W,.^ 0 is the pixel intensity at the reference ROI before 
the flash. The normalized image contrast is also called the 
normalized irradiance or pixel intensity contrast. 

The normalized surface temperature contrast is defined as 
expressed in equations (23), (24), and (25), 


AT-AT„ f 
A T 4- A7V e / 

(23) 

A T re f = T re f - Tfgf, and 

(24) 

AT = T -T°, 

(25) 


where C l is the normalized flash surface temperature contrast, 
AT is the change in the surface temperature of the measure- 
ment ROI after the flash, T is the surface temperature at the 
measurement ROI at the post-flash time t, T° is the surface 
temperature at the measurement ROI before the flash, AT^is 
the change in the surface temperature of the reference ROI 
after the flash, T re/ -is the surface temperature at the reference 
ROI at the post-flash time t, and T re ° is the surface tempera- 
ture at the reference ROI before the flash. 

The normalized image contrast is related to the normalized 
temperature contrast. As a first approximation, an exemplary 
embodiment of the methods described herein assumes that 
the air transmissivity is unity due to short distance (e.g. 1 ft) 
and the air heat radiation is negligible between the object and 
the camera. Therefore, equation (13a) simplifies to equation 
(26), 

W«WVKl-e)WV (26) 

Substituting equations (17) and (18) into equation (26) results 
in equation (27), 

W«C ram o(eT 4 +(l-e)T re /). (27) 

Assuming that the contribution to the pixel intensity from the 
reflected temperature remains the same between the measure- 
ment ROI and the reference ROI, the temperature rise at the 


12 

measurement (object) ROT is related to the rise in the pixel 
intensity at the corresponding image ROI by equations (28) 
and (29), 

AWaC OT „a(€(T 4 -T° 4 )+(l-e)(T^, 4 -T„/ 4 )), and (28) 

AW^C_a(€Cr„/-T„/V(l- 6 )(T„/-T^» 4 )). (29) 

Substituting the equations (28) and (29) in the definition of 
image contrast, equation (20), the method arrives at the 
10 expression in equation (30), 

-g (30) 

( £ ( ( r 4 - 7*) + (7jf - 7$)) + 2(1 - s)(T%f - 7^)) ' 

15 

If a simulation program, such as the software program sold 
commercially known as ThermoCalc, is used to generate the 
temperature evolution at the measurement and reference 
20 ROIs, then equation (30) can be used to generate the image 
contrast evolution provided the method estimates the reflec- 
tion temperature evolution. In practice, the reflection tem- 
perature changes with time due to cooling of the flash lamps. 
The reflected temperature evolution can be estimated based 
25 on the measured data on a high reflectivity foil and a high 
emissivity tape. With the estimation of the reflected tempera- 
ture evolution, the expression of equation (30) can be used to 
estimate the normalized image contrast. 

The estimated image contrast would be partially compen- 
30 sated for the reflection temperature. Only the camera mea- 
surement for the reflection temperature is compensated. The 
reflection temperature evolution is a measure of the heat 
source (or incident) temperature evolution. The net heat trans- 
fer due to the radiation is relatively small. Heat convection is 
35 accounted for in the simulation performed by the Ther- 
moCalc software. The simulation program provides choices 
for the shape and duration of the heat pulse. Depending upon 
the user’s choice, the heat source simulation may or may not 
model the afterglow or the reflection temperature. Thus, if the 
40 reflection temperature is not modeled in the simulation, then 
the estimation of the image contrast from the simulated tem- 
perature would he considered to he partially compensated for 
the afterglow. 

Alternatively, a second approach is to increase the convec- 
45 tive heat transfer coefficient by the average radiation heat 
transfer coefficient. The convective heat transfer coefficient is 
of the order of 5 W/m 2 K for flash thermography. Although, 
the radiative heat transfer coefficient is a function of the 
temperature, it is less than 5 W/nrK for a 10 K difference in 
50 the temperature range of about 300 K ambient temperature. 
The radiative heat transfer coefficient can be lumped together 
with the convective heat transfer coefficient to a value of 10 
W/m 2 K in the simulation. This relationship may allow better 
prediction of the surface temperature in the flash thennogra- 
55 phy using the model and also better prediction of the image 
contrast. 

Yet another approach is to model the heat source as a high 
intensity pulse followed by a very low intensity decaying 
afterglow. 

60 The estimation of the temperature contrast from the IR 
measurement requires measurement of the pre-flash surface 
temperature and the measurement of pixel intensity of a black 
tape affixed to the surface of the test object. For instance, the 
method may use a high emissivity (e>0.95) tape (e.g., tape 
65 sold commercially by the 3M Company under the product 
name of 3M 33 or 3M Scotch Brand 235) positioned near the 
test object. The emissivity of the tape is assumed to be con- 
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stant for both the pre-flash and the post-flash times. The 
method measures the pre-flash surface temperature of an ROI 
in the center of the black tape using a calibrated radiometric 
IR camera. The surface teniperamre needs to be close to the 
equilibrium or steady state. The method also measures the 5 
ROI intensity and assumes a high transmissivity value for the 
air. Thus the pixel intensity due to the tape can be expressed as 
given in equation (31), 

W U,p.=Cc<,m°((l-tt« P e) T re/+£ l <,pe T l <,pJ')- (31) 10 

Due to the high emissivity of the tape, equation (31) can be 
simplified to equation (32), 

^ ' tape’ = ^-' canfi^tape^ 'tape )■ (32) 

15 

If the method measures the pre-flash tape temperature, then 
the method can estimate the camera constant (C' cam ) as given 
in the expression of equation (33), 


tape 


r 7"° 

ape 1 tape 


20 

(33) 


Later a more precise expression for the camera constant will , - 
be derived. 

To measure the T fl or W^j, the method introduces a dif- 
fused high reflectivity (e<0.05 ) reflector such as an aluminum 
foil that is adequately wrinkled. The emissivity of the foil is 
assumed to be constant for both pre-flash and post-flash 30 
times. The foil is placed level with the object and close to the 
physical location of the measurement and reference ROIs, 
such as on the top surface of the test object surface. The 
method can use two foils, one on either side of the measure- 
ment ROI. Here, the method assumes that the reflection tern- 35 
perature is uniform over the test object and the foil, as given 
by equation (34), 

W /oi rC ram o((l -e /oi; )T„/+e*>.7 T f J). (34) 

The contribution of the foil emission is relatively small com- 40 
pared to the foil reflection. In the exemplary method, from 
equation (34), an expression for the reflection temperature 
may be written as expressed in equation (35), 


T 0 4 _ 
4 reft = 



cuu -£/„,,) 


(37b) 


The difference between the fourth power of the reflection 
temperature before and after the flash is now given by equa- 
tion (38), 


t— I pO 4 

‘ reft ~ ‘reft = 


Wjoil-WfoU 

C’caJi-Sjou)' 


(38) 


Knowing the camera constant, equation (3 8) can be evaluated 
as a function of the post-flash time and then substituted in the 
image contrast expression of equation (30). The combined 
steps of using the high reflectivity foil and of using the high 
emissivity tape is referred to as the “Foil-Tape Method” for 
the exemplary embodiments described herein. 

The surface temperature evolution data can be obtained 
using a simulation program such as the previously mentioned 
ThermoCalc software. The IRFT data acquisition using a 
radiometric camera that accounts for the transient reflection 
temperature can provide the surface temperature evolution 
data. Currently, this feature is not available in IR cameras. 
Here, the exemplary embodiments provide an IR measure- 
ment method so that measurement of the transient surface 
temperature can be performed by the IR radiometric cameras. 
The change in the pixel intensity at the foil ROI may be 
defined as expressed in equation (39), 

iWftfW^Wja 0 . (39) 

Using equations (30), (38), and (39), the estimated image 
contrast expression using the highly reflective foil is given by 
equation (40), 


cL = ■ 


£ ((r 4 -ro 4 )-( 


<<t 4 - r" 4 ) + (t^ - O) + 2 

'-'cow ' 1 t'foil) 


(40) 



cut 1 -ejbu) 


45 

(35) A simpler expression for the normalized image contrast will 
now be derived. From equation (13a), the change in the mea- 
surement ROI intensity may be derived as expressed in equa- 
tion (41), 

50 


Also, the exemplary method assumes that the foil temperature 
is about same as the pre-flash reference temperature due to the 
contact with the test object. The black tape is in contact with 
the test object, therefore the temperature estimations of equa- 
tion (36) hold 55 




e(T - t°)(t + r°)(r 2 + r° 2 ) + 

(1 + 1%) 


(41) 


T 1 O^t 0 ^,’T’ 0 

foil= 1 foil = 1 ref tape ■ 

Equation (35) then simplifies to equation (37a), 


(36) 


Similarly, the change in the intensity for the reference ROI 
may be expressed as in equation (42), 



CUT — 44 foil ) 


*( Tret - TrefWref + T^t^ef + T^) + 

(1 -eXT^f, - 1 %XT rtfl + 7^X7^ + T» 2 ) 


(42) 


65 The difference between the change in the measurement ROI 
The pre-flash reflection temperature may then be expressed intensity and the change in the reference ROI intensity may be 

by equation (37b), expressed as in equation (43), 
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e(t- t°)(t + t°)(t 2 + r 02 ) - 

(43) 

r* ~ . 
= 

u 

r (i ) 1 

s(T„ fl - ifefHTref + T%){T 2 rt f + 7^) 

5 


2 - -1 A7^ 

1 + V ’ 

(AT + AT ref ) 


The sum of the change in the measurement ROI intensity and 
the change in the reference ROI intensity may be expressed as 
in equation (44), 


Next, equation (50) may be rearranged for the normalized 
temperature contrast expressed in equation (51), 


AW +AW,ef = 


(44) 


s(T - T°)(T + T t, )(T 2 + T° 2 ) + 


15 


<*J«r - V'fXTr', + tt'fiTl'f + 7 ^) + 
2(1 -eXr«f - T%)(T reJI + + 7^) 


Assuming that the difference between the reference and mea- 
surement temperature is small (e.g. <5 K) at the time of 
maximum contrast, then the expression in equation (45) may 
be shown, 


( t+T°)(t 2 +t 02 ) (4; ” 

= Wnf + + itf) 30 

= (Tnjl + )• 


Therefore, using equations (43) and (44) the expression in 35 
equation (46) may be shown, 


AW-AW ref ^ [e(T-I t> )-s(T ref -T° rtf )} (46) 

AW + AW reJ ~ [e(T - T») + e(T ref - I% { ) + 2(1 - E )(T reJl - T % fl )] ' 40 


Using equation (20), equation (46) may be written as 
expressed in equation (47), 

45 

, _ e(AT - AT ref ) (47) 

w ~ [e(AT + AT r< f) + 2(1 - s)A T rlfl ] ' 

50 

where the change in reflection temperature is expressed in 
equation (48), 

AT^=(T^-T^°). (48) 

55 

Therefore, equation (47) may be written as expressed in equa- 
tion (49), 

(AT-AT ref ) (49) 60 

[(AT + Ar^) + 2^-ljA7'^ ! ] 

Using the definition of the temperature contrast of equation 65 
(23), the method arrives at the normalized image contrast 
expression of equation (50), 




1 + 


2 (;~ 1 K 

(AT + AT r ef) 


cV 


(51) 


Thus, the not finalized image (pixel intensity) contrast and the 
normalized temperature contrast can be approximately 
related to each other. An expression for the normalized pixel 
intensity contrast to normalized temperature (or pixel inten- 
sity/temperature) contrast ratio e' may be expressed as given 
in equation (52), 




A T re ft 


(AT + ATref) 


(52) 


and an expression for the relationship between the normal- 
ized image contrast and the normalized temperature contrast 
may be written as expressed in equation (53), 

CffAeU. (53) 

The pixel/temperature, contrast ratio indicates that the reflec- 
tion temperature rise (caused by flash decay), the object tem- 
perature rise (due to the flash power during rise time), and the 
emissivity affect the conversion between the two contrasts. 
The reflection temperature is due to the glass shield and the 
flash lamps. An estimated amount of less than 50% of the 
flash energy comes through the glass shield of the flash lamps 
during the flash duration. The remaining energy comes 
through the glass shield as the flash lamps continue to glow 
after the flash, providing an increase in the reflection tem- 
perature. Similarly, the glass casing of the bulb and the glass 
shield get hot during the flash and continue to provide thermal 
radiation. Some flash-hood designs use fans to cool the lamps 
after the flash to reduce the amount of heat from afterglow. 
Assuming that the reflection temperature is the same as the 
ambient temperature immediately after the flash, then equa- 
tion (53) may be simplified to the expression in equation (54), 

(54) 

A shutter that covers the flash lamps immediately after the 
flash can be used to achieve a constant reflection temperature. 
In practice, AT rey? is positive and the image contrast and the 
temperature contrast would differ unless the object has an 
emissivity value of one. Now, equations to compute tempera- 
tures at the measurement ROI and the reference ROI using the 
measured ROI intensities will be derived. From equations 
(27) and (31), the temperature at the measurement ROI is 
given by equation (55a), 
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T = 


E tapeT [llpe (1 E )T re fl 


(55a) 


In equation (55a), all quantities are known or measured 
except the reflection temperature. The reflection temperature 
is derived from equation (35) and may be expressed as in 
equation (55b): 


approach to determine emissivity from the IRFT data taken 
with the foil and the tape on the test object starts with the 
relationship between the pixel intensity and the temperature 
for the foil and the tape. Using equation (36), the method in 
5 one embodiment can represent the pre-flash pixel intensity of 
the foil and the tape as in equation (59) and equation (60), 
respectively, 

W yi /<:_o((l-e /o , 7 )T„/+e /o , 7 T^/)aiid (59) 

W^°-C ra „a(e,^T^/+(l-e^)T„/). (60) 



The pre-flash temperature at the measurement ROI is given by 
equation (56), 


Equations (59) and (60) may then be arranged as expressed in 
1 _ equations (61) and (62), 


20 


£_ , jO 4 

.) +I * 


and 


( 61 ) 


^ TV 4 -(!-*</ 

” tape 


(56) 


25 


w° 

W foit 

( 1 - Sfoil ) 



„ 7-e 4 

E foil 1 tape 

(1 — Sfoil ) 


(62) 


Equations (61) and (62) may then be combined, to express 
equations (63), (64), and (65), 


In equation (56), all quantities are known or measured except 30 
the pre-flash reflection temperature, which is given by equa- 
tion (37b). Similarly, the temperature at the reference ROI is 
given by equation (57), 


Tref = 


e^^T-Tl' 4 -( \-s)T re f 

’’tetDe 


35 

(57) 


40 


yv tape "foil 

( 1 — £ t ape ) ( 1 — £ foil ) 




jO 4 rpO 4 ' 

•e 1 tape foil 1 tape 

~ £ tape) (1 ~ £ foit) 


W tape Wfoit _ ^ J £ tape £ foil 4 A 

(1 - e^) ( 1 - s foil ) ~ cam \ (1 - s^) (1 - s foil ) ’ ** 


C T° 4 - r crT 0 4 

'-'cam 1 tape { ~-cam u 1 tape 


w° 

yv tape foil 

( 1 — £ tape ) ( 1 — £ foil ) 

£ tape ^foit 

(1 ~ £ tape) (1 — £ foit ) , 


(63) 


(64) 

(65) 


Similarly, the pre-flash temperature at the reference ROI is 
given by equation (58), 


"t tO 3 


(1 -e)1%a 


45 

(58) 


50 


Equation (65) provides a more precise expression for the 
camera constant defined in equation (33). Equations (59) and 
(60) may then be arranged as expressed in equations (66) and 
(67), 






(1-e,^) ' 


( 66 ) 


Using equation (55a) through equation (58) and equation (23) 
through equation (25), the temperature contrast may be com- 
puted. Thus, the temperature contrast may be calculated using 
the ROI intensity provided by the IR camera. The formula 55 
assumes that the half-max width at the measurement and 
reference point is at least four times the width of the respec- 
tive ROIs so that the ROI. intensity provided by the camera is 
close to the maximum intensity possible for the given tem- 
perature. 60 

The accuracy of the estimation of the temperature contrast 
from the measured pixel intensity data is quite sensitive to the 
value of the surface emissivity of the test object. In one 
embodiment, the method assumes that the emissivity of the 
measurement ROI and the reference ROI are the same. In 65 
another embodiment, the method can either independently 
measure emissivity or use a recommended value. Another 


and 


w%,, c (l )r°/ 

£ foil ^foit 



(67) 


Equations (66) and (67) may then be combined and expressed 
as in equation (68), 


K* „ ((1-^)7-°/ ( 1 - £foii (68) 

= ■ 

£ tape ^ foil £ tape £ foil 


Equations (64 ) and (68 ) may then be combined and expressed 
as in equations (69) and (70), 
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w ° 
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Next, equation (70) may be rearranged and a quantity J may 
expressed as in equation (71), 
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Further, a quantity B may be denoted as expressed in equation , 
(72), 
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Equations (65) and (72) may then be combined and expressed 
as in equation (73), 


(73) 

The camera constant is given by the representation in equa- 
tion (74), 
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The pre-flash temperature at the measurement ROI is given by 
equation (75), 50 

W^C^oteToVl-etT,.,/ 4 ). (75) 

Due to some noted observations and assumptions, including 
contact of the tape with the test object, the steady state of 
surface temperature before flash, and the reflection tempera- 55 
ture not differing more than 5 K from the test object tempera- 
ture before flash, it may be presumed that equation (76) holds 

T°“T,^,°. (76) 

Equation (76) may then be substituted into equation (75), 60 
resulting in the expressions of equations (77), (78), (79), and 
(80), 


W°^C caml r( S Ti; + 


(77) 65 


W" s C a ,„crT" 


20 

-continued 


£ 4 - (1 -£) 


r 0 4 > 
‘reft 


W° =«(s + (l -e)J), and 


( 78 ) 


(79) 


£ ' 



1-7 ' 


(80) 


Thus, the emissivity of the test object at the room temperature 
can be estimated by an embodiment described herein. If the 
measurement is done at two different ambient temperatures 
that are at least 5 K apart, the method may provide better 
confidence in the average emissivity measurement. 

This example illustrates the steps for converting the nor- 
malized image contrast from the thermography data to the 
normalized temperature contrast using additional measure- 
ments required for determining the reflection temperature. In 
one embodiment, the test object may comprise back-drilled 
elongated flat bottom slots or holes at various depths. FIG. 3 
shows the IRFT image of a test object that is made of rein- 
forced carbon-carbon (RCC) with back drilled flat bottom 
slots and with foil and tape (post-flash time=0.6 sec). The 
method in this example is interested in a measurement on the 
center slot. The locations of the measurement and reference 
ROIs are indicated in FIG. 3. One piece of the wrinkled 
aluminum foil is placed on the left side of the test object and 
the other piece of an identical foil is placed on the right side of 
the test object. The locations of the foil ROT are indicated. A 
high emissivity adhesive tape is affixed at the top of the test 
object. The tape ROI is indicated in FIG. 3. The image frame 
time is extremely close to the peak contrast time. The image 
indicates normalized pixel intensity. The image display of the 
contrast is adjusted by manipulating the display rectangle 
location and size. In this example, the measurement ROI is a 
single pixel (or simply lxl pixel). The sizes of the ROIs are 
provided in the Table I. A lxl pixel in dimensions is equal to 
0.023 in.x0.023 in. Normally, ROI cursors and lines have 
colors for quick identification, but colors are not necessary in 
this description. 


TABLE I 


Pixel Size of Region 

of Interests 

Region of Interest 

Pixel Size 

Measurement 

1 x 1 

Reference 

3x3 

Foil 2 

15 x 15 

Foil 1 

15 x 15 

Tape 

9x9 


The dashed cross lines of FIG. 3 are used to plot the pixel 
intensity along the lines to locate the peak point and the 
baseline intensity for the selection of the measurement and 
reference ROIs. FIG. 4 shows a plot of pixel intensity along 
the horizontal pixel line taken in a frame at post-flash time of 
0.6 seconds. Various ROIs, such as the foil, measurement, and 
reference ROI, are also identified in FIG. 4. 

Table II gives the pre-flash measurements and emissivity 
inputs. Note that the representation of Eq. (80) is applied to 
the Table II data (except for the test object emissivity). The 
test object emissivity in this example is estimated to be 0.75 at 
the measurement ROI. A significant number of such measure- 
ments would improve confidence in the value. Therefore, the 
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method may instead use an emissivity value of 0.78, which is 
the value published in the literature for the object material. 


TABLE II 


Pre-flash Measurements and 
Emissivity Inputs 

T„/(K) 

300 

T,^,°(K) 

296 

tape 

0.99 

tfoil 

0.05 

e 

0.78 

W , 0 (average bit value) 

6382 

Wy- oj7 0 (average bit value) 

6768 

W re f (average bit value) 

6480 


FIG. 5 provides plots of the pixel intensity evolutions of the 
two foil ROIs. The vertical dotted line at time equal to 0.6 
seconds indicates the peak contrast time. The foils are reflec- 
tive but specular reflectivity is undesirable. Moreover, the test 
object is also diffusely reflective. The two foils provide two 
independent measurements of the pixel intensity evolution. 
The two intensity evolutions agree with each other and an 
average value is chosen. The two foils are on either side of the 
measurement/reference ROI. Thus, an average of the two foil 
intensity evolutions provides an estimate of foil intensity 
evolution in the measurement/reference ROI with better con- 
fidence as opposed to choosing a single foil intensity evolu- 
tion. 

FIG. 6 shows the pixel intensity evolution at the tape ROI. 
The flash duration is set to 0.003 seconds. The time between 
the image frames is 0.017 seconds for a frame rate of 60 
frames/second. The intensity decays quite rapidly before the 
capture of the first post-flash frame. The test object reaches 
the maximum temperature at the end of the flash duration, but 
the camera is incapable of capturing the corresponding peak 
intensity due to the slow frame rate. FIG. 6 indicates that the 
temperature rise is higher on the tape ROI than that on the foil 
ROI. The tape ROI is in the area of tape overhang and indi- 
cates similar intensity evolution as that of the foil ROI. 

FIG. 7 shows the pixel intensity evolution at the measure- 
ment (upper curve) and the reference (lower curve) ROI .FIG. 
7 indicates the hot spot (higher intensity) of the measurement 
ROI in comparison to the reference ROI. The two traces for 
the measurement ROI and the reference ROI eventually 
merge together as the temperature of the indication area and 
the surrounding area become equal. 

FIG. 8 shows a plot of a simple image contrast, which is 
defined as the difference between the intensities of measure- 
ment and reference ROIs . The simple contrast evolution at the 
measurement ROI indicates a peak time of approximately 0.5 
seconds. 

Referring now to FIG. 9, this embodiment begins with 
equation (65) being used to compute the camera constant. 
Next, equations (37a) and (37b) are used to compute the 
reflection temperature evolution. FIG. 9 provides an example 
of a computed reflection temperature evolution. The pre- flash 
reflection temperature is estimated to be about 305 K, which 
is 5 K above the pre-flash temperature of the test object. Also, 
the reflection temperature is 2 K higher than the pre-flash 
reflection temperature. When the data for this example was 
taken, the setup was used to acquire five other shots within 5 
minutes before this represented shot, and the time between 
the last shot and die current one was less than a minute. Thus, 
the flash lamps were warmer than the ambient temperature 
before the flash and caused the reflection temperature to be 
higher than the ambient temperature. In at least one of the 
embodiments described herein, the method computes the 


measurement temperature evolution using equation (55a). 
Using equation (57), the method also computes the reference 
temperature evolution. The three estimated temperature evo- 
lutions are shown together in FIG. 10. The upper curve is for 
5 the reflection temperature. The middle curve is for the mea- 
surement ROI temperature and the lower curve is for the 
reference ROI temperature, FIG. 11(a) shows the image con- 
trast computed using the pixel intensity and one fourth (1/4*) 
power of pixel intensity to illustrate equation (45), which 
10 suggests that the contrast computed with the fourth power of 
the temperature and the contrast computed with unity power 
of the temperature are approximately the same. FIG. 11(b) 
demonstrates the estimated temperature contrast (upper 
15 curve) along with the image contrast (lower curve) with 
€=0.78. Thus, the pixel and temperature contrast differ, if the 
test object’s emissivity is less than one. FIG. 11(b) also illus- 
trates the difference in the peak contrasts due to lower values 
of the emissivity. FIG. 12 shows an exemplary screen shot of 
20 the steps of computing the temperature contrast evolution 
(upper curve) and pixel (image) contrast evolution (lower 
curve) with software implementing the IR Contrast method 
which relates the thermal measurements in terms of the nor- 
malized anomaly contrast as a function of the frame number. 
25 FIG. 13 represents the plot of the contrast ratio e'. Two plots 
are shown. One plot (lower curve) uses the representation 
given by equation (52). The other plot is the ratio of the two 
contrasts shown in FIG. 12. The plots differ slightly at early 
instances of time, but the two values agree well at later times. 
30 The contrast ratio is approximately the same as the emissivity 
of the test object for longer times when the rate of change in 
the surface temperature is relatively low. 

In at least one embodiment, the method estimates the nor- 
malized image contrast based on results from the surface 
35 temperature simulation of the software sold commercially as 
ThermoCalc in order to understand the effect of the reflection 
temperature on the profile of the pixel contrast. In this par- 
ticular example, the method takes a sample simulation tem- 
perature evolution at the measurement and reference areas. 
40 The simulation power input is 1 .8E6 W/nT2with a duration 
setting at 0.003 second. This set-up assumes about 30% effi- 
ciency for converting the electrical power (12 kJ source) to 
the absorbed energy. Depending upon the flash power setting, 
the reflection temperature evolution is affected. The method 
45 of this example estimates the reflection temperature evolution 
based on the experimental method explained in this descrip- 
tion. 

FIG. 14 shows the simulated temperature evolutions with 
an assumed reflection temperature evolution. The upper curve 
50 represents the reflection temperature. The middle curve is 
indicative of the measurement ROI temperature and the lower 
curve is for the reference ROI temperature. The normalized 
temperature contrast is calculated using the representation 
given by equation (51). The method of this particular example 
55 assumes that the emissivity of the test object equals 0.78. 
Using equation (30), the method of this example can compute 
the image contrast. The results are shown in FIG. 15 . The 
upper curve shows the temperature contrast. The lower curve 
shows the image contrast. Using equation (27), the pixel 
60 intensity can be determined if the camera constant is known. 
Thus, this example illustrates the influence of the emissivity 
change on the pixel contrast. It also suggests that changes in 
the pulse power have a small influence on the pixel contrast 
due to changes in the reflection temperature. 

65 If the method of this particular embodiment assumes that 
one can input the ambient reflection temperature in the simu- 
lation model, then the afterglow of the heat source would be 
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due to differences in the radiation between the pre-flash and 
post-flash reflection temperature. The afterglow flux may be 
represented by equation (81), 

(«) 5 

The afterglow flux evolution starts from the first post- flash 
frame. The flash lamp pulse flux evolution is estimated for a 
duration between the start of the flash and the first post-flash 
frame. An example of afterglow flux evolution is given in 
FIG. 16. FIG. 17 shows an example of flash flux evolution to 
until the first post-flash frame. The flash lamp pulse and the 
post-flash afterglow evolutions can he added to model the 
compound source flux evolution for the entire time of data 
acquisition. 

FIG. 18 is a flow diagram illustrating methods for convert- 15 
ing an image contrast evolution of an object to a temperature 
contrast evolution, for modeling afterglow flux, and for 
assessing the test object emissivity, in accordance with some 
embodiments described herein. Processing begins at step 
1800 whereupon, at block 1805 , a steady state temperature 20 
measurement is made of the tape affixed to the test object 
before the flash or the application of heat. At block 1810 , the 
transient thermography IR digital video data of the object, 
foil, and tape is collected. The collection of data begins before 
the application of heat. At block 1815 , the pixel intensity 25 
evolutions at the measurement, reference, tape, and foil ROIs 
are measured. At block 1820 , the normalized pixel intensity 
(image) contrast evolution is computed using equation (20). 

At block 1825 , the camera constant is computed using equa- 
tion (65) and the reflection temperature evolution is computed 30 
using equation (37). At block 1830 , the measurement tem- 
perature evolution is computed using equation (55a) and the 
reference temperature evolution is computed using equation 
(57). At block 1835 , the normalized temperature contrast 
evolution computed using equation (23) or equation (51 ). At 35 
block 1840 , the afterglow flux is modeled using equation 
(81). At block 1845 , the test object emissivity is assessed 
using equation (80). Processing subsequently ends at step 
1899 . 

FIG. 19 is a flow diagram illustrating a method for convert- 40 
ing a temperature contrast evolution of an object to an image 
contrast evolution, in accordance with some embodiments. 
Processing begins at step 1900 whereupon either, at block 
1905 , the source flux and the reflection temperature evolution 
are assumed, or at block 1910 , the source flux is based on the 45 
heat source setting and modeled from the experimental 
reflection temperature evolution using equation (81). At 
block 1915 , the transient thermography surface temperature 
on a void-like flaw in a selected material is simulated. At 
block 1920 , the measurement ROI temperature evolution and 50 
the reference ROI temperature evolution are simulated. At 
block 1925 , the normalized temperature contrast evolution is 
computed using equation (23), followed by block 1930 where 
the normalized pixel intensity contrast evolution is computed 
using equation (30) or equation (50). Alternatively, at block 55 
1935 , the camera constant is calculated from equation (33) 
and the change in pixel intensity evolution at the measure- 
ment ROI and the reference ROI are computed using equation 
(28) and equation (29), respectively, followed by block 1940 
where the normalized pixel intensity contrast evolution is 60 
computed using equation (20). Processing subsequently ends 
at step 1999 . 

Hie normalized image contrast and the normalized tem- 
perature contrast differ for objects with emissivity other than 
one. Therefore, for accurate results, the embodiments and 65 
examples described herein indicate that the two quantities 
should not be treated as the same. To compare the simulation 
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temperature contrast with the measured pixel contrast, the 
method of the embodiments described herein should estimate 
the evolution of the reflection temperature and the incident 
heat flux. Ideally, the set of instructions from a computer 
software program that implements the methods described 
herein should model the compound heat source flux evolu- 
tion, which also includes the thermal afterglow. The effect of 
the reflection temperature on the pixel intensity should also be 
accounted for to seek a better estimation of the temperature 
contrast evolution from the pixel intensity evolution data. 

As described herein, the reflection temperature evolution is 
established based on data acquired by the IRFT data acquisi- 
tion system. This system acquires readings from a test object, 
a high emissivity tape with known emissivity, and a diffuse 
highly reflective metal foil with known reflectivity. The 
method also records the steady state pre-flash temperature of 
the object using a thermocouple (or other contact sensor) or 
an accurate radiometer. The IR datacube is recorded using the 
normal IRFT technique. 

The method comprises four regions of interests (ROIs). 
One region is for the measurement ROI. The second region is 
for the reference ROI. The third region is the foil ROI for 
measurement of the reflection temperature. The fourth region 
is for measurement of the pre- flash temperature at the high 
emissivity tape. 

Using formulas described previously herein, the method 
estimates the reflection temperature evolution. The method 
then computes the temperature contrast from the IRFT data. 
In accordance with at least one embodiment described herein, 
the pixel intensity/temperature contrast ratio has been defined 
in order to relate the temperature contrast to the image con- 
trast. 

The embodiments described herein also include a method 
for using the evolution of the reflection temperature to model 
the afterglow flux of the flash source. Using the estimated 
compound source evolution in simulation software, the tem- 
perature contrast evolution may be estimated and then the 
image contrast profiles on the simulated voids may be esti- 
mated. 

The embodiments described herein include an emissivity 
estimation technique using the same IR camera that is used in 
an IRFT system. This method provides determination of the 
emissivity for the desired thermal wavelength. By using the 
foil-tape technique during the IRFT shot, the transient reflec- 
tion temperature or the reflection temperature evolution can 
be recorded. If an IR camera is programmed with the repre- 
sentative formulas for reflection temperature formulas 
described herein, the camera can provide the object surface 
temperature directly, even during the IRFT data acquisition. 
The IR camera can also be programmed to estimate the object 
emissivity in real-time by using the formulas described herein 
in combination with the set-up described for the foil-tape 
technique. 

In light of the principles and exemplary embodiments 
described and illustrated herein, it will be recognized that the 
exemplary embodiments can be modified in arrangement and 
detail without departing from such principles. Also, the fore- 
going discussion has focused on particular embodiments, but 
other configurations are contemplated. In particular, even 
though expressions such as “in one embodiment,” “in another 
embodiment,” or the like are used herein, these phrases are 
meant to generally reference embodiment possibilities, and 
are not intended to limit the invention to particular embodi- 
ment configurations. As used herein, these terms may refer- 
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ence the same or different embodiments that are combinable 
into other embodiments. 

Similarly, although exemplary processes have been 
described with regard to particular operations performed in a 
particular sequence, numerous modifications could be 5 
applied to those processes to derive numerous alternative 
embodiments of those described herein. For example, alter- 
native embodiments may include processes that use fewer 
than all of the disclosed operations, processes that use addi- 
tional operations, and processes in which the individual 10 
operations disclosed herein are combined, subdivided, rear- 
ranged, or otherwise altered. 

In view of the wide variety of useful permutations that may 
be readily derived from the example embodiments described f . 
herein, this detailed description is intended to be illustrative 
only, and should not be taken as limiting the scope of the 
invention. What is claimed as the invention, therefore, are all 
implementations that come within the scope of the following 
claims, and all equivalents to such implementations. In the 
claims, means-plus-function and step-plus- function clauses 
are intended to cover the structures or acts described herein as 
performing the recited function and not only structural 
equivalents, but also equivalent structures. Thus, while a nail 
and a screw may not be structural equivalents in that a nail .,,- 
employs a cylindrical surface to secure wooden parts 
together, whereas a screw employs a helical surface, in the 
environment of fastening wooden parts, a nail and a screw 
may be equivalent structures. 

30 

The invention claimed is: 

1. An apparatus for converting an image contrast evolution 
of an object to a temperature contrast evolution, the apparatus 
comprising: 

one or more processors; and 35 

one or more memory units coupled to the processors, the 
apparatus being configured to: 

calculate a measurement region of interest temperature 
change AT; 

calculate a reference region of interest temperature 40 
change AT,. e/ ; 

calculate a reflection temperature change AT,. e7? ; 
calculate the image contrast evolution U,,, 7 ; and 
convert the image contrast evolution to the temperature 
contrast evolution C 1 according to the equation of: 45 


c' ~ 


1 + 



(AT + A7V <f ) 


K 


50 


wherein 6 is the emissivity of the object. 

2. The apparatus of claim 1 further comprising: 

an IR camera operatively connected to the one or more 55 
processors. 

3. The apparatus of claim 2 wherein the apparatus is further 
configured to measure and display the object surface tem- 
perature in real-time. 

4. The apparatus of claim 2 further comprising: 60 

data acquisition electronics coupled to the IR camera; 

a flash lamp; 

a power supply/trigger unit coupled to the flash lamp; and 

a computer interfaced with the power supply/trigger unit 
and the data acquisition electronics, the computer being 65 
configured to trigger the flash lamp and initiate data 
acquisition. 


5. The apparatus of claim 4 further comprising: 

a test object positioned at the focus of the IR camera; 
a fell attached to the surface of the test target; and 
a tape attached to the surface of the test taiget. 

6 . An apparatus for converting an image contrast evolution 
of an object to a temperature contrast evolution, the apparatus 
comprising: 

one or more processors; and 

one or more memory units coupled to the processors, 
the apparatus being configured to: 

calculate a measurement region of interest temperature 
change AT; 

calculate a reference region of interest temperature 
change AT,. e/ ; 

calculate a reflection temperature change AT,. e/? ; 
calculate the image contrast evolution U, r 7 ; and 
convert the image contrast evolution to the temperature 
contrast evolution where the apparatus being config- 
ured to calculate the image contrast evolution com- 
prises: 

measuring a measurement region of interest tempera- 
ture T; 

measuring a reference region of interest temperature 

Ue/i 

measuring a reflection temperature T re /, and 
calculating the image contrast evolution according to 
the equation of: 


c* w 




' (£((7-4 - 7U) + (7^ - T? ef 4 )) + 2(1 - sXT?^ - 7°/)) 


wherein e is the emissivity of the object, t° is the measure- 
ment region of interest temperature at time of flash, T re ° 
is the reference region of interest temperature at time of 
flash, and T re;/f 0 is the reflection temperature at time of 
flash. 

7. An apparatus for converting an image contrast evolution 
of an object to a temperature contrast evolution, the apparatus 
comprising: 

one or more processors; and 

one or more memory units coupled to the processors; 
the apparatus being configured to: 

calculate a measurement region of interest temperature 
change AT; 

calculate a reference region of interest temperature 
change AT re/ ; 

calculate a reflection temperature change AT^; 
calculate the image contrast evolution C^; and 
convert the image contrast evolution to the temperature 
contrast evolution, where the apparatus being config- 
ured to calculate the image contrast evolution com- 
prises components configured and arranged for: 
measuring a measurement region of interest tempera- 
ture T; 

measuring a reference region of interest temperature 

T refi 

measuring a reflection region of interest pixel inten- 
sity W fol7 ; 

calculating a camera constant C' com ; and 
calculating the image contrast according to the equa- 
tion of: 
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wherein 6 is the emissivity of the object, € rbil is the emis- 
sivity of a foil, T° is the measurement region of interest 
temperature at time of flash, T re f is the reference region 
of interest temperature at time of flash, and W^,, 0 is the 10 
reflection region of interest pixel intensity at time of 
flash. 

8. The apparatus of claim 7, where the camera constant is 
calculated according to the equation of: 
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wherein s. tape is the emissivity of the tape, T tape ° is the tape 
region of interest temperature at time of flash, and W tape ° 
is the tape region of interest pixel intensity at time of 23 
flash. 

9. An apparatus for converting a temperature contrast evo- 
lution of an object to an image contrast evolution, the appa- 
ratus comprising: 

one or more processors; and 30 

one or more memory units coupled to the processors, the 
apparatus being configured to: 

calculate a measurement region of interest temperature 
change AT; 

calculate a reference region of interest temperature 35 
change AT,. e/ ; 

calculate a reflection temperature change AT,.^; 
calculate the temperature contrast evolution C ; and 
convert the temperature contrast evolution to the image 
contrast evolution, 40 

where the apparatus being configured to convert the tem- 
perature contrast evolution to the image contrast evolu- 
tion is calculated according to the equation of: 



(AT 4- A T ref ) 


wherein 6 is the emissivity of the object. 

10. The apparatus of claim 9 further comprising: 

an IR camera operatively connected to the one or more 55 
processors. 

11. The apparatus of claim 10 wherein the apparatus is 
further configured to measure and display the object surface 
temperature in real-time. 

12. The apparatus of claim 10 further comprising: 60 

data acquisition electronics coupled to the IR camera; 

a flash lamp; 

a power supply/trigger unit coupled to the flash lamp; and 

a computer interfaced with the power supply/trigger unit 
and the data acquisition electronics, the computer being 65 
configured to trigger the flash lamp and initiate data 
acquisition. 
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13. The apparatus of claim 12 further comprising: 

a foil for attachment to the surface of the object; and 
a tape of known emissivity for attachment to the surface of 
the object. 

14. An apparatus for converting a temperature contrast 
evolution of an object to an image contrast evolution, the 
apparatus comprising: 

one or more processors; and 

one or more memory units coupled to the processors, 
the apparatus being configured to: 

calculate a measurement region of interest temperature 
change AT; 

calculate a reference region of interest temperature 
change AT,. e/ ; 

calculate a reflection temperature change ATj^; 
calculate the temperature contrast evolution C; and 
convert the temperature contrast evolution to the image 
contrast evolution, where the apparatus being config- 
ured to calculate the temperature contrast evolution 
comprises components configured and arranged for: 
measuring a measurement region of interest pixel 
intensity W; 

measuring a reference region of interest pixel inten- 
sity 'Nrefi 

measuring a tape region of interest pixel intensity 
W 

” tape j 

measuring a reflection temperature T,.^; and 
calculating the temperature contrast evolution 
according to equation of: 


AT^AT^ 

AT+A7V 

a r= r-r° 


wherein 


A T„, = T, 


ref ~ ■* ref 1 ref ■ 


Kf, T = 
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r . 4 „ 4 
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"ref Q 4 _ . j. 4 
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and T Pe j = 


p re $ T° 4 _ 1 1 _ p1T° 4 

c-tape y-Q 1 tape ' A refl 
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wherein 6 is the emissivity of the object, e tape is the emis- 
sivity of the tape, T tape ° is the tape region of interest 
temperature at time of flash, iff’ is the reflection tem- 
perature at time of flash, W,.^ 0 is the reference region of 
interest pixel intensity at time of flash, and W tape ° is the 
tape region of interest pixel intensity at time of flash. 



